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Abstract:  A general  emission  equation,  which  contains  as 
its  asymptotic  limits  both  the  Richardson  (thermal 
emission)  and  the  Fowler  Nordheim  (field  emission) 
equations,  is  developed.  The  general  thermal-field 
equation  performs  over  a range  of  temperatures,  fields,  and 
work function  values  where  the  asymptotes  are  degraded. 
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Introduction 

Formulae  for  the  evaluation  electron  emission  current 
density  are  generally  predicated  on  the  high-field  low- 
temperature  (Fowler  Nordheim  Equation)  or  low-field 
high-temperature  (Richardson-Laue-Dushmam  Equation) 
limits  of  the  general  emission  integral  given  by 


m,Pr)  = ^ \D{E,)f{E,)dE, 

/(£,)=  ^?^ln[l  + exp(^r(/^-^.)]  [1] 
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large.  For  all  conditions  of  interest,  the  large  u limit 
suffices,  and  is  dropped  from  further  consideration.  The 
integral  in  Eq.  [2]  can  be  separated  into  regions  such  that 
series  expansions  of  the  components  of  the  integrand  can 
be  integrated  term-by-term  exactly.  Collecting  terms  in  the 
resulting  series  expansions  gives  rise  to  the  leading  order 
approximations  to  N(n,s)  of 
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The  regions  n < 1 and  n > 1 shall  be  referred  to  as  the 
“thermal”  and  “field”  regimes,  respectively.  Eq.  [2]  using 
Eq.  [3]  constitutes  the  Generalized  Thermal-Field  model. 
The  “revised”  FN  and  RED  equations  become 

JpN^  ^RLD  ^ iPp - ■E’o ))Z {Pf ! Pt'\ 

J RLD  '^RLD  {PbPt  ) ^ {Pt {M~  EgifL  {Pf  / Pp  } 


where  D(E)  is  the  linearized-0  form  of  a hyperbolic-tangent 
approximation  to  the  transmission  probability,  f(E)  is  the 
supply  function,  fip  and  fir  = are  slope  factors,  F 
[eV/nm]  is  the  product  of  the  field  and  the  electron  charge 
e,  p [eV]  is  the  chemical  potential  (Fermi  level),  y[A/cm^] 
is  current  density  and  other  terms  have  their  usual 
meanings.  E^  is  the  energy  related  to  the  component  of 
electron  momentum  directed  at  the  surface  - henceforth, 
the  subscript  shall  be  suppressed.  Eg  is  reference  point  to 
be  determined,  and  the  coefficient  C is  of  order  unity  (and 
will  be  taken  henceforth  as  unity).  Eq.  [1]  can  be  recast  as 


where  the  inherent  symmetry  between  the  two  limits  is 
manifest  and  the  argument  of  E deserves  particular  notice. 
Note,  however,  that  when  n is  approximately  unity,  Eq.  [4] 
is  no  longer  adequate  and  the  n « 1 approximation  must  be 
used.  The  remaining  tasks  are  to  determine  E„  and  fp. 

For  a potential  barrier  described  by  V(x),  the  WKB 
approximation  to  the  “area  under  the  curve”  0is  given  by 

0{E)=  \2  (V(x)  - Ef^  dx  [5] 


J{Pf,Pt)  = ApldT^N  {E,  - p\PpE, 
-Pp 

N(n,s,u^=  n j Inj^l  + dz 


[2] 


where  the  limits  of  the  integral  are  at  the  zeros  of  the 
integrand  and  V(x)  > E otherwise.  If  we  define  E„  to  be  the 
location  of  the  integrand  maximum  in  Eq.  [1],  then 
expressions  for  E„  and  ftp  are  given  by 


where  = 120.17  A/cm^K^.  The  ratio  n = is  of  Pp  ~ = E^y,  Eg  = E,,, +(d(E^y  fip^  [6] 

central  importance:  its  behavior  governs  the  transition  from 
thermal  to  field  emission  as  n goes  from  very  small  to  very 
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For  E > V(x),  a linear  extension  of  0(E)  is  taken:  from  the 
parabolic  potential  WKB  approximation  we  find 

where  Q = 0.36  eV-nm  and  (/>  = <i>-  ^40^  . In  searching 
for  the  integral  maximum,  is  an  elusive  quarry  that  is 
best  hunted  for  numerically.  For  thermal  emission 
conditions,  it  is  found  near  the  barrier  maximum,  whereas 
for  field  emission  conditions,  it  is  found  near  the  chemical 
potential.  For  Cu-like  parameters  (//  = 7.0  eV,  <?=  4.6  eV, 
T=  800K),  its  behavior  is  shown  in  Figure  [1]. 


Figure  1 Location  of  current  integrand  maximum  for 
copper  parameters  at  800  Kelvin. 


Figure  2 Behavior  of  the  parameter  n as  a function  of 
field  for  copper  parameters  and  T = 800  K. 

Moreover,  it  has  been  found  that  the  parameter  n satisfies  a 
power-law  relation  with  respect  to  F at  fixed  T.  Using  the 


parameters  of  Figure  [1],  Figure  [2]  shows  the  behavior  of 
n cc  E^.  The  value  of  the  power  p is  obtainable  through 
simple  models,  which  shall  be  taken  up  separately. 


Figure  3 Comparison  of  Eqs.  [3]  & [4]  to  commonly 
used  forms  of  the  FN  and  RLD  equations  for  copper 
parameters  but  for  a low  O of  2.4  eV  and  T = 1000  K. 

The  performance  of  the  Generalized  Thermal-Field 
equation  (GTFE)  is  shown  in  comparison  to  standard 
implementations  of  the  Fowler  Nordheim  and  Richardson 
equations  but  for  “challenging”  parameters:  low  work 

functions  at  moderate  temperatures.  The  GTFE  results  are 
in  agreement  with  numerical  evaluation  of  Eq.  [1].  Apart 
from  the  expected  disagreement  in  the  n = 1 region,  it  is 
also  seen  that  standard  FN  result  falls  at  high  fields  for  low 
work  functions,  and  that  standard  RLD  overestimates  at 
low  fields  because  of  the  impact  of  £’„’s  shifting  location. 

In  summary,  a generalized  Thermal-Field  emission 
equation  has  been  derived  based  on  a linear- 
approximation  for  the  potential  evaluated  at  the  maximum 
of  the  current  integrand.  The  transition  from  “thermal-like” 
emission  to  “field-like”  emission  is  demarcated  by  the 
transition  of  n = PEPf  from  <1  to  >1.  from  near  the 
barrier  maximum  to  near  the  Fermi  level.  The  method  has 
general  utility  for  potentials.  Revised  FN  and  RLD 
approximations  were  given  and  shown  to  give  good 
agreement  with  the  exact  method. 
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